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Abstract 

The superconducting properties of the 2D fermion system with local and different types of the 
indirect boson-exchange attractions in the cases of s-wave and d-wave pairing are reviewed and 
analysed at T = 0. In particular, the possibility of the crossover from the Bose-Einstein conden- 
sation regime to the BCS-like superconductivity with the carrier density and coupling changing in 
the case of different pairing channels is discussed. Gaussian fluctuations of the order parameter 
are taken into account and the carrier density dependence of the gap is studied in this case. The 
role of the interaction form in the physical behavior of the systems is also discussed. 
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I. INTRODUCTION 



The theoretical description of the carrier density dependence of the physicsl properties of 
superconductors has a long history. Probably, the first attempt to solve this problem self- 
consistently was made by Eagles in Ref. where the author tried to apply his results for 
description of the superconducting properties of Zr-doped SrTiOs P|. The author studied 
the dependence of the superconducting gap at T = and of the mean-field critical temper- 
ature on charge carrier density in the two-dimensional (2D) and in the three-dimensional 
(3D) systems with phonon-like attraction at low carrier densities. The set of the equations 
for the superconducting gap and Fermi energy in two and three dimensions was analyzed. 
It was estimated that the diameter of pairs is smaller than distance between pairs at low 
carrier densities, therefore, superconductivity at low carrier densities transforms into super- 
fluidity of small (local) pairs. Now this phenomena is known as crossover from superfluidity 
to superconductivity with carrier density changing. It was also shown by Eagles that there 
exists a critical value of the inter-particle coupling below which there is no local pairs at low 
carrier densities in the 3D system. It is worth to mention that the possibility of superfluidity 
scenario of superconductivity was also proposed by Ogg 0| and Schafroth 0,0,0 before 
paper Q]. ^ 

Later the problem of the crossover was studied by Leggett in |t7||, where the author 
studied 3D system with short-range repulsion and finite-range attraction. He analyzed the 
properties of the system changing dimensionless parameter l/lakp), where a is the scattering 
length and kp is the Fermi momentum which defines the particle number in the system. It 
was shown that in the limiting cases the system consists of bound "diatomic molecules" at 
l/^akp) = oo, and the Cooper pairing takes place at l/^akp) = — oo. 

Nozieres and Schmitt-Rink generalized the results of Leggett on the cases of finite 
temperatures and lattice model when the interaction potential is separable. They have 
shown that the crossover from superfluidity to superconductivity is smooth with kp and 
interaction (attraction) potential V changing, when V is larger than the critical value for 
two-particle bound state formation. All the results mentioned above correspond to the case 
of the s-wave pairing. 

The real boom in the interest to the crossover phenomena has started in the second half 
of eighties after the discovery of high-temperature superconductors (HTSCs), materials with 
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very unusual dependence of superconducting properties on the carrier density. 

Theory of the HTSC still remains one of the most difficult and one of the most important 
problems of the modern condensed matter theory. Due to complicated crystal structure, 
anysotropic pairing, strong electron correlations, presence of the disorder, low dimensionality 
of this materials and some other complications, the construction of the self-consistent theory 
of HTSCs is not completed so far. 

During the last years many models which take into account some of the cuprate properties 
were proposed for describing the doping and interaction dependence of the superconducting 
gap and of the superconducting critical temperature. 

The 3D crossover from superfluidity to supercondictivit y in the g-wave pairing channel 
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1611 for the model 



for the model with local attraction was studied in 

m 

with local (so called 4F-) attraction, and in for the case of separable potential. In ^ 



'he role 



(2D 



the additional case with local repulsion and inter-site attraction was considered, 
of the order parameter fluctuations at T = was analyzed in [18] (3D case) and in 
and 3D cases). The quasi-2D model with s-wave pairing at T = was considered in [21 
For the 2D case the problem of the crossover in the s-wave pairing channel was analyzed 



at zero temperature (w^ 




26 



U> 15, 17, 



len a long-r ang e superconducting order is still possible in a 2D system 
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and at finite temperatures in 
29|. The metal-superconductor phases boundaries on the n — U phase 
diagram of the 2D attractive Hubbard model at T = was studied in by means of the 



Dynamical Mean-Field Theory. The crossover problem for the 2D phonon-mediated model 
in the s-wave channel at T=0 was studied in [31! [3^. 

The 2D c/-wave crossover for models with separable potential was studied in 
at T = 0. The Ginzburg-Landau potential with carrier density dependent parameters for 
the 2D (i-wave superconductor at finite T was considered in 3^. The doping dependence of 
the critical temperature in the strongly correlated electron model with the electron-phonon 
coupling was studied in It was shown that the vertex correction to the electron- 

phonon interaction in this model leads to the increasing of superconductivity in the d-wave 
pairing channel. In paper |2^ the cases with different pairing symmetries s, s^y, d^y and 
dx-2-yy in the 2D system with local, nearest neighbor, and next nearest neighbor attraction 
were studied at T = and finite T in order to describe the doping dependence of the 
zero superconducting gap and superconducting critical temperature. The possibility of zero 
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temperature s-, d- and mixed s + id pairing in the 2D system as function of coupling were 



studied in paper [37[. The same problem in the quasi-2D Hubbard model with nearest 
neighbor attraction at finite temperature was studied in [s^ , where also doping dependence 
of the superconducting properties was considered. 

The 2D model with the nearest neighbor (n.iO attraction and next nearest neighbor 
(n.n.n.) hopping at T = was considered in 40|. It was shown that at some relation 
between n.n. and n.n.n. hopping the system is always in Bose-Einstein condensation regime 
and that there is no pairing at low carrier density when coupling is weak (see also |33[). 
However, as it was shown in |^| this statement is not correct. The 2D model with local 
repulsion and n.n. attraction at T = in the s-wave and rf-wave channels was studied in 42 1. 



The role of the n.n.n. hopping was analyzed. The s-wave and d-wave pairing crossover at 
T = in the model with doping dependent attraction was considered in 43[. It is interesting 
to mention that the boson-fermion model with electron and hole sectors and different kinds 
of fermion-boson coupling was proposed in ^ to unify the BEC and BCS phenomena. 
Superconductivity or superfluidity regimes take place at different model parameters. 

Despite the d-wave pairing symmetry is considered as a typical property of cuprates 
4^, there are experimental evidences that the s-wave or mixed s + d pairings take place 
in some materials at certain values of doping. Thus, it was shown that optimally doped 
Pri s^^Ceo,i45Cu04^_y at low temperatures demonstrates nodeless gap inconsistent with d- 
wave pairing and electron doped superconductor SrnQLan iCuO-?. demonstrate s-wave 



pairing near optimal doping 



43] • The analysis by Zhao j4S| shows that dominant bulk sym- 



metry of the order parameter in some cuprates is extended or anysotropic s-wave symmetry. 
It was also found that the crossover from rf-wave pairing in underdoped and optimally doped 
regime to s+ci- wave pairing in overdoped regime takes place in YBa2Cu307-s [49]. Crossover 
from (i-wave pairing to s-wave pairing with doping near optimal c was found for electron 
doped Pr2-xCexCuOi [50j and for Pr2-xCexCu04^-y, La2-xCexCuOi-y |51j. However, the 
s-wave pairing is not always present in overdoped cuprates, d-wave regime in overdoped 
Tl2Ba2Cu30e+s was found in p^ . 

Below we analyze the behavior of the superconducting systems with different inter-particle 
potentials in the s-wave and d-wave pairing channels at T = as functions of particle density 
and interaction strength. We shortly review the main properties of different systems. We 
apologize for these authors whose results are not presented here, since it is very difficult to 
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cover all the results in such wide and fastly growing field as theory of superconductivity, 
even in a case of some special topic. 



II. THE MAIN EQUATIONS 

The most general Hamiltonians, which is usually studied in theory of superconductivity 
can be written as: 

H — — y ^ tnmCna^^ma ~ A* ^ ] ^na^na ~ ^ ^ ] "^It "^l J, "^n J, Cnf — ^ ] Vnm(^n}^^mi^'^i^n-\ 
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n,q n,q 



1) 



_m{q) 

where Cna = c^ir, n) is a fermionic field operator with spin a I at lattice site n and at 
time r, tij describes the nearest neighbor (n.n.), next nearest neighbor (n.n.n.) and other 
hopping processes; fi is the chemical potential of the system. The inter-particle interaction 
is described by the terms proportional to Vq (local interaction) and Vnm (n.n. or n.n.n. 
interaction). The last two terms describe additional electron boson interaction and free 

— * 

boson parts of the Hamiltonian, where g is a boson mode with coordinate Xn{q), momentum 
Pn{q) and frequency cUniq) and gn{q) is the fermion-boson coupling. One can easily pass to 
the continuum version of this Hamiltonian with substitution of the n.n. hopping operator tnm 
by tnm Snmt{l — (a^/(2(i)) V^) (where a is the inter-cite distance, d is the dimensionality 
of the system), introducing radius cut off in the the interaction terms, etc. 

In the case of the d-dimensional square lattice the free fermion dispersion relation in the 
momentum space has the following form, when the n.n. hopping takes place 

d 

C{k) ^ -2tJ2cos{akj) - 1^, (2) 

j=i 

where k is a d-dimensional wave vector. 

It is convenient to find the thermodynamic potential by using the path integral approach 
for studying the properties of a quantum many-particle system. This method is not necessary 
in the case of the mean-field solution, but it is extremely useful when the fiuctuations are 
studied. 

The partition function of the system is 

Dip^Dipe-^ (3) 



with the action 



S 



dr 



(4) 



To study the superconducting properties of the system, we make the Hubbard- 
Stratonovich transformation with bilocal fields (pnmi'Tii ^2) and ''"2) can be apphed 
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exp 



Km(Tl, T2) 

-0L(n,r2)^„i(ri)V^„l(r2) -7/'i^(ri)^i^^(r2)0„™(ri,r2)) , (5) 

where V^m(Ti,T2) includes the effective inter-particle attraction due to boson couphng after 
intergations over the field X. 

Let us introduce the Nambu spinor 

^n(r) = f "^f^^ I ,^1{t) = «(r),^„,(r)) . 

In this case the partition function can be written as 

where 

= [ dri [ dT2 y^{M^!!iililIl)J ^(^^^ _ ^^^^t (^^^) - f^(t„m - ^nmyU)] ^m(l"2) 

Jo Jo ^ Vnm{Ti,T2) 

-0L(n, T2)^i{n)f_^mir2) - ^Un)T+^„(r2)0„„(ri, r2)}(7) 

where f± = ^{f^ ± fy) are the Pauli matrices. 

The last action is diagonal on fermionic fields, therefore the integration over \E'''" and 
can be performed exactly. The partition function becomes in this case 
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where fl[Q] is the thermodynamic potential, which in the "leading order" is 

pn[g] = r dn r dr, y '^"^^'""'1'' - TrLnG-^ + TrLnG^\ (8) 

Jo Jo t;;^ VnUn,r2) 

The Nambu spinor Green function G satisfies the following equation 

= 5{Ti- T2)6nm 

with anti-periodic boundary conditions 

Gnmin - T2+ f3) = -Gnm{n - T2) . 

The thermodynamic potential (jH)) is the most general form of the superconducting effective 
potential with non-local retarded inter-particle interaction. It will be used below to study 
the fluctuation effects. 

The minimization of the thermodynamic potential with respect to the order parameter 
and chemical potential leads to the following system of the coupled equations for ^nmiTn, Tm) 
and fi: 

=0 (9) 



or 

$nm(l"l, T-2) = Vrd{r\, T3,)trf+Glm{r3, T2) (11) 



n/ = -trf^G„„,(ri,r+), (12) 

where Uf = Nj/v is particle density in the system {v is the volume of the system). 

In general it is very difficult to find the Greens function Gnm(Ti,T2), therefore some 
simplifications must be applied. We shall consider the case of the space and time invariance, 
Vnm{Ti, T2) = Vn-m,o{Ti — ^2) ■ In this casc the Green's function has the following form in the 
momentum space: 

iujn + r3^(k) - k)f+ - k)f_ 



G{iun, k) 
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and the system of the equations (fTTj). (fT^ has the following form: 



V{p,k)+gl 



(2^ 



E 



$(t^m, P) 



UJt. 



e(p) 



|'^'(t^m,P) 



^(p-k)^ 



(27r)^-(^™-a;„)2-a;(p-k)2j 



(13) 



(27r)" 



e(p) 



e(p)2- |<l>(^^,p) 



(14) 



V(p — k) is the Fourier transform of the non- retarded interaction and the term proportional 
to Qb describes the inter-partice boson coupling, uJn = TiT{2n + 1) is the standard Matsubara 
frequency. Interaction term in (jl3|) is written in general form, it describes, for example local 
non-retarded interaction, when \/(p — k) = const, = 0, non-local non-retarded interaction, 
when V{p — k) 7^ const, gb = 0, local retarded interaction, when V{p — k) = const, gj, ^ 
0,uj{p — k) = const etc. The system fll3|) . ()14|) will be analyzed in the next Section for 
different forms of the inter-particle potential V(p,k) and spectra ci;(p — k). 

III. SOLUTIONS 

A. Model with local non-retarded attraction 

The problem of the crossover from small to large fermion density in the model with 
ocal attraction was considered in for the 3D case and in 
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27| for the 2D case. 



For the simplest case of local non-retarded attraction the interaction parameters in equa- 
tions (Uni), (fT^ have the following form: V(p, k) = V = const, gb = 0. Therefore, the gap 
in this case is momentum and frequency independent 



$(a;„, k) = A = const. 

The summation over frequency in (fT^ . (fT^ can be easily performed, and one gets the 
standard system of the equations: 



1 = V 



tanh 



/V^(k)2 + A2^ 



(27r)'^2v/C(k)2 + A2 \^ 2T 



(15) 
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(2^ 



e(k) 



tanh 



which have a simple form at T = 0: 



V 



(2vr)'^2v/e(k)2 + A2 



(27r)« 



1 



e(k) 



(16) 



(17) 



(18) 



Since the gap is momentum independent, the only isotropic s-wave pairing regime is possible 
in the model with local attraction. 

The s-wave pairing regime in the case of quadratic dispersion law 

k2 



e(k) 



2m 



and indirect pairing with boson energy cut-off 6[f|^(p) — /i| — cUc) was considered in the 2D 
and 3D case for low carrier densities by Eagles In fact, in this case the integration over 
momentum j d'^k/{2TiY can be replaced by integration over energy j p{e)de, where p(e) is 
the density of states (DOS), it is constant in the 2D case and ~ y/e in the 3D case. 

It is very easy to solve system (|T7jl and (|TH|l in the 2D case. In the case when ep << W 
{W is a free fermion bandwidth) the solution has a simple form: 



A ~ v^2W^e-2-/(™^), 



where €b = —2We~^'^^^"^'^^ is the two-particle bound state energy. Obviously, the crossover 
from superfluidity to superconductivity with doping takes place at any coupling constant, 
there exists the value of ep when /i = for any V. 

This is not true in the 3D case, when the crossover takes place only when the coupling 
constant is larger than some critical value Vc- This difference follows from the difference in 
the DOS in the gap equation. In the 2D case p(e) = const and the gap equation has the 
solution A = A/2Vr|£:b| at /i = and any coupling constant. 
In the 3D case the gap equation has the following form 

k'^dk 1 



1 = V 



27r2 2Vk^ + A2 
9 



(19) 



at /i = 0. The integral over k on the right hand side has maximal value \/2mW /(47r^) at 
A = 0. Therefore, a simple estimation for Vc is given by the relation 1 ~ VcV2mW /(47r^), 
or Vc ~ 4:n^/V2mW. 

In the case of the momentum cut-off / (|^6'(|^(p) — /u| — cOc) the approximate solution 
can be also easily obtained in the 2D case: 

A ~ ^/2\eb\ABcsdi^^c - cf) + /^Bcsd^ep - uj^) 

where A^cs = 2uJce~'^'^^'^'^^^ is the BCS expression for the gap and the bound state energy 
in this case is Eh = —2ujce~^'^^^"^^\ The gap is growing with doping in this case, and 
asymptotically reaches its maximal value A = Abcs when ep » Uc- 

In the 3D case the DOS in equations (fTTj) and (fTHj) can be substituted by the DOS on 
the Fermi level, and the 3D solution can be obtained from the 2D one with substitution 
mV/{2iT) kprnV / {Att"^). In is possible to estimate the critical value of the coupling 
constant when the crossover takes place: Vc = 27r^/ y/muJ^. In both cases with and without 
momentum cut off the transition from superfluidity to superconductivity is smooth, in other 
words the gap is groving with doping continuously. 



B. Models with non-local non-retarded interaction 



It is important to study a more realistic case of the non-local attraction in presence 
of the short range Coulomb repulsion. In order to study the superconducting properties 
of such a model in the channels with different pair angular momentum /, it is convenient 
(see 29, Q, 

for example), to approximate the interaction potential by by a separable 

function: 

Vl(ki,k2) = -A;«7KkiVKk2), (20) 
where A; is an effective coupling constant, and 



wi{k) = hi{k)cos{lip]^) 



(21) 



hi{k) 



{k/k,y 

[l + k/koY+'Z^' 



(22) 
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k is the momentum modulus = |k| and y^k is the momentum angle in polar coordinates 
k = k{cos{(p]i), sin{ip]i)). Parameters ko and ki put the momentum range in the proper 
region - the potential is attractive at rg < r < ri and repulsive at r < rg, where ko ~ I/vq 
and ki ~ l/ri. 

It is easy to see that the interaction potential (|2(jp has the correct asymptotic behavior 
at small and large momenta: Vi(kik2) ~ k{k2 and Vi(kik2) ~ 1 / V^i ^2 , correspondingly. 
Since the region of low carrier concentrations, where the crossover can take place, is the 
most interesting, the correct behavior of the interaction potential at small momenta is the 
most important. The small momenta give the main contribution to the integrals in the case 
of low carrier concentrations (see equations below). We shall study s- and (i-wave channels 
with / = and 2 separately, so we assume that the parameters A; for both channels are 
independent. 

In this case the equations for the gap and for the chemical potential have the following 
form: 

dp Ai{p) 



A^k) = -A; 



dp 



i2nr2^e\p)+Mp} 



-Mp,k), 



(2vr) 



'e'ip) 



The solution of equation has the following form 

AKk) = Af^^w^ik), 



(23) 



(24) 



(25) 



where A[°'* does not depend on the momentum k. 



As it was shown in 
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29|, the crossover from superfluidity to superconductivity is 



smooth with doping increasing (see Figs. 1 and 2 below). However, there exists a critical 
value of the interaction potential in the (i-wave pairing channel, below which the crossover 
from superconductivity to superfluidity is impossible in the (i-wave pairing channel. 



The more realistic case in connection to HTSCs was considered in |43|, where the corre- 
lation length To was considered at small carrier densities as tq ~ a/ ^Jnj. This dependence 
for the length of spin-spin correlations at small carrier densities was found for example in 
La2-xSrxCu04^. The magnetic correlation length decreases with carrier density per cell in 
this material as 3.8A/ ^/nj [54]. The value of a was considered to be a = A/2/7rao, oq is the 



square lattice constant. The following equality gives the value for a: {ji /2)rlNf = alNceii, 
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where on the left side the volume of the 2D system is expressed as a volume (circle of the 
radius ~ ro) occupied by one particle, multiplied by the full number of particles Nf, Nceii is 
an elementary cell number in the system. The free fermion bandwidth W is connected with 
ttQ as W = TT'^/{mal). It should be noted, that the relation rg ~ a 7,^/^7 at a = ^J2pKaQ is 
in a good agreement with the experimental data for La2-xSrxCuOi [^, where the plane 
magnetically ordered lattice parameters are equal to 5.354 A and 5.401 A, and the corre- 
sponding parameter a is ^ 3.8 A. 

It was shown the the critical value of the coupling exists even in the s-wave pairing 
channel for this case (see Figs.l and 2). 
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FIG. 1: Crossover line coupling-carrier density is presented for the s-pairing channel (solid line). 
The dotted line represents the corresponding curve for the case ro(nj) = const at rg = oq. The 
insert shows the doping dependence of the crossover value for coupling at very low charge carrier 
densities Here and below all parameters are expressed in units of the bandwidth W 



Different versions of the previous model can be considered. Namely, the correlation radius 



ro{nf) = a/2/ (vrn^) can be introduced in a model with exponential decay of the attraction: 



Vb{r,t)=g, 

Its Fourier transform has the following form: 



-r/ro 



e-b- 



A(q) 



ge-bTo 



(27r)Vl + roV' 
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FIG. 2: The crossover line coupling-carrier density is presented for the d-wave case (solid line). 
The dotted line is the crossover curve for the case ro(nj) = const at ro = ao 



In addition, the same kind of the short-range repulsion can be introduced in the model. It 
is easy to see, that in the 2D case this potential has the following form (see, for example 
0): 

^e(q) = ^ ^7,,, y (28) 

where g^-e = 27re^, the Thomas-Fermi momentum is qtf = 4e^m/7r = A/TcaB, and qb = 
l/(e^m) ~ 0.529 Ais the Bohr radius. The function g{x) is defined as 
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(x) = i-e{x- 1)^/1- 1/ 



(29) 



This model in the s-wave pairing channel demonstrates the crossover from superfluidity 
to superconductivity at any value of the coupling constant, contrary to the previous case, see 
Fig. 3. The doping dependence of the gap and of the chemical potential at different values 
of the dimensionless interaction parameters Xe-b = gl-b^'^ol -^e-e = 5'e-e"^^o/(47r) 
is presented in Figs. 4 and 5. As it follows from these Figures, the gap is decreasing with 
doping growing. This situation is in a qualitative agreement with the experiments for the 
cuprates. One can consider also an interesting situation of the combined local-|-non-local 
attraction, when the local attraction will tend the Cooper pairs to transform into the local 
pairs. 

The "mixed" case with hi{k) = 1 and A = + Adcos(20k) was considered in Ref. 
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FIG. 3: The crossover interaction value A^™^* (defined by the condition fi = 0) is presented as a 
function of ej? at Xe-e = 0. 
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FIG. 4: The gap (a) and the chemical potential (b) as functions of ep at different Xe-b and 
A._. = 0. 



IstI Issj l . In particular, it was shown in that the crossover from the d-wave to the s-wave 
superconductivity with the intermediate s + d-phase ta kes p lace with doping in agreement 
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with the experiments on some cuprate materials ^ , ^ 

The d -wave case W|^(k) — cos{kx) — cos{ky) was considered in 

HHQ- In particular, 
the role of the next nearest neighbor hopping t' on the on the pairing was studied in 0, 0| . 
As it was shown in |41|, the crossover with doping take place at any value of t', when the 
coupling constant is bigger than V^- 
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FIG. 5: The gap (a) and the chemical potential (b) as functions of ep at different Xe-e and 

Xe-b = 0.01. 

It was stated in [2^, that the model with small on-site repulsion and nearest neighbor Vi 
and next nearest neighbor V2 attraction with V2 ~ 60 — 80meV and V2/V1 ~ 1.3 — 1.5 can 
describe the experimental data for hole-doped oxides. 



C. Models with retarded interaction 



The retardation effects in the interaction can play very important role the supercon- 
ducting properties of the system. Let us consider a boson propagator with the dispersion 
u = u;(k): 

uj^ — uj^[k) + id 

where the spectra a;(k). For example, in the general case of phonon dispersion one has 
a;(k) = a/coIq + c^k^. We shall study this rather general case below with the following 
approximation: D{u!, k) ~ D{u!, k^?), so the effective propagator can be written as 

where >/uj^'+~c^k^ is substituted by new effective frequency Wq- This approximation corre- 
sponds to the case of the optical phonon attraction (in the case when ujq ^ 0). Generally 
speaking, the gap is frequency-dependent in this model. The problem of the crossover in the 
model with frequency dependent gap was studied in [s^ . 

The system of the equations for the gap and for the chemical potential in this case is 

<P[U) = -ig J (27r)2 _ ^(k)2 _ |0(j,)|2 + ,s{uj- z/)2 - ooi + l5' ^ ' 
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2eF = Re / duj[{— 

'0 



+ 1M^-V/^2_|^(^)|2)]. 



(33) 



Let us describe shortly how this system of the equations can be analyzed analytically. 
First of all, it is possible to show that the approximation <^{uS) = A = const in PHj) is a 
rather good one. Then, this equation gives: 

\e,\ 

where the two-particle bound state energy e\j depends on the coupling parameter in this case 
(see below). 

After the Wick rotation uj iuj the gap equation becomes 

d^kdv 6(iy) 1 



(uj) = g luq 



(27r)3 z/2 + (k2 /2m - /i)2 + (j)^{u) {lo - vf + ul ' 
In polar coordinates, after integration over the angle we obtain: 

kdk 



(34) 



0(^) 



„2, ,2 /•+00 

9 ^0 



du- 



(27r)2 {u - vf + ujI Jo (P/2m - /x)^ + z/2 + 02(j,) " 

Since 0(c<j) is an even function of u, it can depend only on cu^, we restrict the integration 
over u to the positive values: 



dv(j){v) 



TV a 
— h arctan — — 

2 V^^2 + ^2(^) 



2 Jo Vz/2 + 02(zy) (a;-i/)2 + ^5 
The coupling constant A = g'^m/ {27[) is introduced. The asymptotic behavior for ^(0;) is 



(35) 



0(t^) luj^const-* 0, (f){uj) luj^oo^ — . 

tiai0,0: 



We make the next replacement in the interaction potential 



1 



1 



{uj - z/)2 + a;2 ^0-2 + 
Then, the differentiation with respect to u gives: 



ducp^u) 



(z/2 + ^2)2 + 02(^) 

It is evident that 0' < 0, i.e. (pmax = 0(0) = A. 



/i 



1 H — arctan — , 

TT Vz/2 + 02(^) 



(36) 



(37) 



(38) 
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Differentiating once more and introducing a new variable x = z^^/i^q S^^^ following 
differential equation 

2 A 



x+1 



b'ix) + 



4^{x + iy^/x + {(f){x)/ujoy 



1 H — arctan — -z 

rf ^x + (0(x)/cjo)2 



0. 



(39) 



0, 



(40) 



with the boundary conditions 

<P'{x)l^, = 0; [<^(a;) + (a: + l)</.'(^)]U 

These conditions follow directly from (jSSI), and (jHEI)- 

The analysis of this equation shows that the approximate solution for the gap in both 
weak and strong coupling regimes is (see js^] for details) 



AM 



UJ 



where A is a parameter, which depends on the coupling and the carrier density, what is very 
important. In the weak coupling regime the parameter A has the following coupling and 
carrier density dependence: 



A = ^/2\eb\uJo6{uJo - ep) + AscsOiep - ujq) 
where Abcs = 2a;oe~^/'^ is the BCS expression for the gap and the bound state energy in 



this case is Eb = —2ujQe_ '^^^ 
large carrier densities 



31 



32 



n the strong coupling regime — A and /S.bcs ~^ (4/3)Act;o at 



To summarize, the crossover with carrier density and coupling changing in this effective 
model with retarded interaction is smooth and the gap is not small when the pair frequency 
is smaller than the boson frequency. This approach can be used for studying the crossover 
for the cases of different symmetries of the order parameter. 

It is important to mention that the case considered here can be generalized on the case 
when one takes into account the vertex correction to the electron-phonon interaction. This 



correction is usually small when ep » ujq (the Migdal theorem). It was shown in 



58|, that 



this correction is also small when ep << (^o- However, in some cases the vertex correction 

QBQ. 



can lead to strong enhancement of superconductivity 



]. It is also necessary to note. 



that this correction even lead to enhancement of the ci-wave superconductivity in the strongly 
correlated electron system jsSjIs^, despite phonon interaction due to the symmetry does not 
allow (i-wave pairing in the case when the short range electron repulsion is not considered. 
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IV. FLUCTUATIONS 



The role of fluctuations of the order parameter in the 2D and even in the 3D case at 
T = can be significant. As it is shown in the Gaussian fluctuation correction to the 
s-wave order parameter is non-negligible even in the weak-coupling case. It was shown that 
the fluctuations of the order parameter phase lead to increasing of the gap. 

In this Section we show that the simultaneous order parameter modulus and phase fluc- 
tuations in the model with local attraction lead to strong increasing of the order parameter 
at small carrier densities and to small decreasing of the order parameter when the carrier 
densities are large. 

At zero temperature the thermodynamic potential of the system with local attraction has 
the following form: 

\^ r (pk 

In other words, it depends on the combination |$p = (3?$)^ + ($5$)^. For studying fluctu- 



/(0(vmT^-e(k)) 



ations of the OP, it is convenient to use new real variables: 



<Pi{x) 

such that = (j)\{x) + (t)\{x) = |$(a;)p. Another possibility is to choose the decompo- 

sition of the order parameter on its phase and modulus $(a;) = A{x) exp{i6{x)). However, 
it leads to some difficulties, since one needs to keep the order parameter modulus positive 
at functional integration over the fluctuations. The "old" order parameter variables are 
connected with the new ones in the following way: 

$(a;) = (P^(x) + 202(x), = 01 (x) - #2(3;). 

We assume that the mean-field value of the field ^ is chosen as 

A 



%{x) = 

The order parameter $ can be written as 







6(f)2{x) 
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in the case when its fluctuations are considered. We neglect the fluctuations of the carrier 
density ^/(r), and consider the homogeneous constant value of over the lattice nf{r) = 
rif = const. 

Substitution of (j41|) into the expression for the thermodynamic potential gives the fol- 
owing correction to the thermodynamic potential in the second order in the fluctuations 



5VL = - 



dv 



(27r) J (27r)' 



6^{iu, q)A{iu, q)5$(--iz/, -q) 



(42) 



where A(zz/, q) = + xl^'^? q) is a 2 x 2 matrix with the susceptibility components 



-tr 



d^k 



2 J.^{27t)J (27r)2 
where lj± = a; ± z//2, k± = p ± q/2. Integration over u can be easily performed: 



Xi2[ii^, qj 



X22[ii', qj 



d^k 1 E+ + E2 



(27r)2 2i/2 + + E_)2 _ 

d^k 1 



1 + 



-X2iU«^,qj = 
d^k 1 



(27r)2 2 + (E, 



E+ + E2 



(27r)2 2z/2 + (^^ + ^_)2 _ 



E_)^]E+E_ 
A2' 



E^E_ 



where E± = + ^2, e± = (k ± q/2)2/(2m) — /i is the free fermion dispersion relation. 
After the integrating out of the fluctuation fleld 5> the correction to thermodynamic potential 
has the next form 

""^ du f d^q ' 



5n = -, 

2 . Loo (27r) 



(2 



:ln 



7r 



+ Xu{i^, q)}{^ + X22{iJ^, q)} - Xi2(«i^, q)X2i(^z^, q) 



The factor ~ A which appears in the measure of functional integration over A in Z due 
to taking into account the symmetry of the thermodynamic potential with respect to the 
transformation $ 6*°"$ (see, for example [g^) is omitted in the last expression. This 
factor can be absorbed in the measure of the functional integral over A in the partition 
function, where the functional integration can be actually performed over the variable A2. 
Let us note, that only the flrst component {y + Xii(^'^)q)) under the logarithm in 6Q will 
be present when one considers the particular case of the order parameter phase fluctuations. 

It is useful to diagonalize the matrix A, to flnd the contributions which come from the 
phase and the modulus fluctuations of the order parameter. Obviously, the flrst component 
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will correspond to the phase fluctuations and the second - for the modulus fluctuations, 
as it follows from the deflnitions of the fleld $. So, it is easy to arrive to the following 
representation 



6Q 



2J_^(2ti)J (27r) 



■In 



{-^ + Xe{ii^, + Xa(«z^, q)} 



where 



Xe{iv,(\) = ^ [xuiiiy,q) + X22(«z^,q)] - xl i [Xiil^'^.q) - X22(«z^,q)]^ - Xi2(iz^,q)X2i(««^,q) 



1 



xaI^z^, qj 



2 [xu{i^, q) + X22(^'^, ^)] + \l\ [xii{i^, q) - X22(^'^, q)]^ - xMii^, q)x2i(«'^, q) 



are the effective contributions to the thermodynamic potential from the fluctuations of the 
order parameter phase and the modulus, correspondingly. 

The equations for the gap and for the chemical potential and (fTUI) have the next 
form in the case of the Gaussian correction to the thermodynamic potential due to the order 
parameter fluctuations: 

du f d'^q [ 9x0(iz/, q)/9A ^XAi'>'l^,^)/^A 



A 

V 



d'q 



A 



1 

+ 7^ 



(27r)2 2^^(q)2 + A2 2 



n 



f 



d\ 



1 



A2 



{2Ti)j (27r)2 [l/r + xe(^//,q) l/V + xi^{i^,^)\ ' 
°° du r d^q \ dxe{i^,(l)/dfi ^ dxA{i^,c[)/dn 



1/F + Xe(^^^,q) l/V + XAiiiy,q) 



,i2n)J (27r)2 

Considering the limit of weak coupling, the functions xe ^^^^ Xa in the denominators 
can be neglected, than the integration over u can be performed and after some trivial 
manipulations the following simple equation can be obtained: 



1 

V 



dq 



V 



(27r)2 2^e(q) + A2 2 J {27, f J {27,f ElE. 



d'^k 



(43) 



dq 



Ve(q) + 



V 



d'^k 



d^q A2(^+-^_ 



(27r)2 J (27r) 



(44) 



(compare with (|17|) and (|18p). The substitution q — >■ — q in a part of the terms was made 
in obtaining equations (j43|) and (j44j) . It is interesting to note that in the case of the phase 
fluctuations the numerator under integral in the last term of the the gap equation will be 
i^_i_(^_ — ^+) and the last term in the number equation will be multiplied by |. It is also a 
good approximation to put fj, = ep, since in the weak coupling regime /x is different from ep 
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FIG. 6: The dependence of A on ei? is presented for the mean-field solution case (solid line) and for 
the order parameter fluctuations case (dashed line) at different values of the dimensionless coupling 
parameter G = mV/{2Tr). The dotted line is the estimation of Kos and Millis for the case of 
the order parameter phase fluctuations. 



only at very low carried densities. The solution of the equation ()4H|1 for the gap parameter as 
function of ep fi = ep and different values of the coupling constant is presented in Fig. 6. 
The estimation of the order parameter in the case of phase fluctuation is also presented. As 
it was shown in in the 2D case the phase fluctuations lead to an effective increasing of the 
coupling constant V V{1 + The gap can be calculated from the standard mean- field 
BCS equation. The comparison of different cases shows that the phase fluctuations lead to 
the gap increasing, while the total fluctuations lead to much stronger increasing of the gap 
at small carrier densities and to decreasing of the gap when the carrier density is large. The 
last result is familiar, but the first one is very surprising. The dependence of the gap on 
coupling at a small value of the carrier density is presented in Fig. 7. These results suggest 
that the higher order corrections in fluctuations should be studied to better understand the 
behavior of the system. 

It should be mentioned that the role of the disorder due dopants in the fluctuations of 
the inhomogeneous order parameter was studied by Yu.G. Pogorelov and by the authors 



m 



62 



63l l64| . We don't describe this important problem here, since it deserves a special 



review. 
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0.2 



0.15 - 




FIG. 7: The dependence of A on G is presented for mean field solution case (solid line) and for 
order parameter fluctuations case (dashed line) at ep = 0.1. The dotted line is the estimation from 
[l^ for the order parameter phase fluctuations. 

V. CONCLUSIONS 

In this paper the crossover from superfluidity to superconductivity with doping increasing 
at T = in the cases of the s-wave and d-wave pairing was reviewed. In the 3D case this 
crossover does not take place at small couplings, and the same situation takes place in the 
(i-wave pairing case in two-dimensions, when the interaction does not depend on the doping. 
In the case when the correlation radius depends on doping the minimal value of coupling for 
the two-particle bound state exists even in the s-wave channel. Also the gap can decrease 
with the doping in this case. 

It was also shown, that the fluctuations of the order parameter play an important role 
at T = 0. The fluctuations of the order parameter phase in weak coupling limit in the case 
of the s-wave pairing regime lead to increasing of superconductivity at any physical carrier 
density, while the modulus fluctuations lead to much stronger increasing of superconductivity 
at small carrier densities, while at large carrier densities they lead to suppression of the order 
parameter, as a result the gap is decreasing in BCS regime, when both the modulus and 
the phase fluctuations are taken into account. This means that the higher order fluctuation 
corrections must be investigated in order to understand the self-consistent theory at low 
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carrier densities. 

We would like to mention some direction which can be interesting for the future investi- 
gations. Studying of the problem of the crossover with realistic dispersion relations is not 
performed even on the mean-field level in many interesting cases. The interplay between 
disorder and superconductivity, and strong correlation and superconductivity is another 
interesting topic for the investigations. 

The fluctuations in the (i-wave pairing channel and in other non-isotropic pairing channels 
even in the case of Gaussian fluctuations were not studied carefully so far. It is also important 
to go beyond the Gaussian fluctuations, since the pair susceptibility is diverging in the 2D 



and in the 3D cases, as it was mentioned in Ref . [19| . The role of the inter-layer coupling is 
another problem which is not solved in general at the moment. The solution of the problem 
mentioned above will lead to better understanding of the superconducting properties of 
systems with arbitrary carrier density and pairing potential. 
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